More on the narrowing of impact broadened radio recombination 
lines at high principal quantum number 



M.B. Belli 



ABSTRACT 

Recently Alexander and Gulyaev have suggested that the apparent decrease in impact broad- 
ening of radio recombination lines seen at high principal quantum number n may be a product of 
the data reduction process, possibly resulting from the presence of noise on the telescope spectra 
that is not present on the calculated comparison spectra. This is an interesting proposal. How- 
ever, there are serious problems with their analysis that need to be pointed out. Perhaps the 
most important of these is the fact that for principal quantum numbers below n = 200, where 
the widths are not in question, their processed generated profile widths do not fit the widths of 
the processed lines obtained at the telescope. After processing, the halfwidths of the generated 
and telescope profiles must agree below n — 200 if we are to believe that the processed generated 
linewidths above n — 200 are meaningful. Theirs do not. Furthermore, we find that after apply- 
ing the linewidth reduction factors found by Alexander and Gulyaev for their noise added profiles 
to our generated profiles to simulate their noise adding effect, the processed widths we obtain 
still do not come close to explaining the narrowing seen in the telescope lines for n values in the 
range 200 < n < 250. It is concluded that what is needed to solve this mystery is a completely 
new approach using a different observing technique instead of simply a further manipulation of 
the frequency-switched data. 
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The theory of impact broadening of radio 
recom bination lines was developed by Grieml 
(Il967l). A broader coverage w as carried out by 
Gordon and Sorochenkol (2002) and radio recom- 
bination lines have been a valuable tool for many 
years for the purpose of studying the physical 
conditions inside galaxies ( Gordon 2008[) . Using 
the frequency switching observing technique we 
have found evidence that the impact broadened 
linewidths of hydrogen recombination lines near 
6 GHz appear to become much narrower than 
predicted at high quantum numbers, n > 200 
(lBel]||l997tlBeU et al. |[2000l:lBeU et al.ll201ll: iBeU 
20111 ). Although a possible explanation for this 



wa s sugge s ted b y [Oki (l2004l) . this was refuted 
bvlGriemI (l2005l). It was later demonstrated by 



Belli (j2011l ) that when observations at other radio 
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frequencies were taken into account the line nar- 
rowing appeared to be correlated with the density 
of reco mbination li nes in frequency space. More 
recentlv iHevI (2012) has obtained results that may 
eventually provide a theoretical explanation for 
the discrepancies seen between what the impact 
broadening theories predict and what has been 
observed in these high-n Rydberg-Rydberg re- 
combination spectra from Galactic HII regions. 

Because the multiple overlap data reduction 
technique we used requires s pecial processing, it 
has re cently been suggested bv lAlexander and Gulvaev 
(j2012i ) that the narrowing at high n may be re- 
lated to the processing. After generating theo- 
retical impact broadened line profiles for Orion, 
these authors claim for the zone between An = 11 
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to 14 (near 6 GHz) that after noise is added to 
the generated profiles their processed hnes exhibit 
a l ine narrow i ng si milar to what was reported 
by iir et alj as being a mystery. Here 

we examine their treatment of the data and find, 
unfortunately, that there appear to be some dis- 
crepancies that make their conclusion much less 
convincing than it appears to be at first glance. 

2. Analysis 

In Fig 1 we have plotted as filled circles the 
processed linewidths we obtained for the Orion 
lines observed using the NRAO 1 40-ft telescope 
(|Bell et al. I I2OOOI: iBeh et all l2Qll[ ). We used a 
frequency-switched offset of ±4 channels which re- 
sulted in a separation of 8 channels between sig- 
nal and reference lines. In processing the data 
we applied 5 overlaps before cleaning. When fre- 
quency switching with multiple overlap reduction 
is used it is important to realize that the processed 
widths of lines located at the edge of the observ- 
ing window cannot be trusted since they can be 
contaminated by edge effects, especially in strong 
continuum sources like Orion. For this reason we 
have previously not included the Orion (n. An) = 
(194,7) transition in linewidth plots and for this 
reason it has also not been included here in Fig 1. 

In Fig 1 the solid curve represents the pro- 
cessed linewidths obtained for our generated im- 
pact broadened Orion Voigt profiles. For transi- 
tions below n — 200 the processed linewidths ob- 
tained for the generated profiles fit the processed 
linewidths obtained for the telescope lines remark- 
ably well for the (102,1), (129,2), (147,3), (184,6), 
and (202,8) lines, all of which are located in rea- 
sonably unconfused regions of the spectrum. The 
width of the (174,5) line is wider than expected 
and this is assumed to be due to the fact that the 
(286,24) line is located in its high-frequency wing. 

Above n — 200 the widths of the processed tele- 
scope lines quickly become much narrower than 
the widths of the processed generated profiles. 
One systematic effect that might possibly have ex- 
plained this line narrowing was the fact that, un- 
like the telescope profiles, our generated profiles 
did not contain a superimposed noise component. 
However, it seemed very unlikely that this would 
introduce such a rapid narrowing and it was con- 
cluded that the Orion lines must become intrinsi- 




Fig. 1. — (filled circles) Widths obtained for our 
telescope lines plotted vs n. (solid curve)Widths 
obtained for our generated Orion Voigt profiles, 
(dashed curve)Widths of Orion profiles generated 
by Alexander and Gulyaev. 




Fig. 2. — Generated Orion linewidths plotted vs 
71. (solid curve) Widths for our generated Orion 
Voigt profiles, (dashed line) Impact broadened 
linewidths calculated by Alexander and Gulyaev 
for Orion. 
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Fig. 3. — (solid and dashed curves) same as Fig 
1. (dotted curve) Processed widths obtained by 
Alexander and Gulyaev for generated Orion pro- 
files after adding noise. 




Fig. 4. — Processed Orion widths plotted vs n. 
(soHd curve) same as Fig 1. (dot-dash curve) 
Widths of our processed profiles after adding 
noise, obtained using the (noise/no noise) reduc- 
tion factor found by Alexander and Gulyaev. 



cally narrower than predicted above n — 200. 

Also plotted in Fig 1 as a das hed curve are the 

proces sed linewidths obtained bv lAlexander and Gulvaev 
(j2012l ) for their generated Orion profiles be- 
fore adding noise. These two curves should 
be identical. Clearly they are not, with the 
dashed curve increasing much more slowly than 
the solid one. These results can be explained 
if either, a) the Orion pro files generated by 



Alexander and GulvaevI (|2012r ) are different from 



the Voigt profiles we generated or, b) their pro- 
cessing has been carried out differently. But there 
are other problems. Their linewidths, given by 
the dashed curve, are not a good fit to the tele- 
scope data below n = 200 where the widths are 
not in question. In addition to this, their curve 
already shows signs of becoming narrower at high 
n- values, even though there has been no noise 
added to the generated profiles at this point. 

In Fig 2 the widths of our generated Orion Voigt 
profiles before processing are shown by the solid 
curve. They can be comp ared to the before pro- 
cessin g widths generated bv lAlexander and Gulvaev 
( 2012f ) shown by the dashed curve. Clearly here 
as well there is a significant difference between 
the two curves with the widths of their generated 
profiles increasing much more slowly than ours. 
This difference then appears to be the explanation 
why the two curves in Fig 1 do not agree. How- 
ever, since the processed widths of our generated 
curves are in good agreement with the processed 
telescope lines below n = 200 while theirs are not, 
we can assume that it is our impact broadened 
Voigt profiles that have been generated correctly. 

This presents a serious problem. If their pro- 
cessed widths below n = 200 do not fit the widths 
obtained for the processed telescope lines in this 
n-range how can we believe that the widths they 
obtained for n > 200 lines are reliable? Further- 
more, the processed widths of their generated pro- 
files are already starting to narrow at high n even 
though there has been no noise added to the spec- 
trum. This indicates that at least some of the 
narrowing they find at high n is unrelated to the 
presence of noise on the spectrum and may be due 
either to the shape of their generated profiles or 
their reduction process. 

In Fig 3 th e dotted curve represents the w idths 
obtained by [Alexander and Gulvaev ( 2012( ) for 
their processed profiles after adding noise to their 
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generated spectra. The solid and dashed curves 
are the same as in Fig 1. Above n = 220 the 
noise added widths begin to decrease rapidly and 
it was concluded by these authors that this ex- 
plains the decrease seen in the Orion telescope 
data plotted here as filled circles. This might be 
a valid conclusion but it would have been much 
more convincing had their processed generated 
linewidths agreed more closely with the processed 
telescope linewidths below n — 200. 

Our calculated Orion Voigt profiles were gen- 
erated by L. W. Avery almost 15 years ago. Un- 
fortunately these profiles, along with our reduc- 
tion programs, are no longer unavailable to us. 
Otherwise it would have been easiest simply to 
add noise to those profiles and reprocess them. 
In an attempt to obtain an approximate noise- 
added result for our processed generated profiles, 
in Fig 4 we have applied the same re d uction fac- 
tor found bv [Alexander and Gulvaev (2012), be- 
tween their before and after noise-added curves 
in Fig 3, to our processed generated widths. The 
result is shown by the dot-dashed curve in Fig 
4. It assumes that all of the narrowing is due 
to added noise, which may not be the case. The 
line narrowing can be seen to begin at significantly 
higher n than does the narrowing in the telescope 
linewidths and it does not agree well with the 
widths obtained for the telescope data, especially 
in the range 200 < n < 240 (8 < An < 14) 
where Alexander and Gulyaev agree that noise 
does not overwhelm the line profiles obtained at 
the telescope. In fact, in the n- region where a 
fit is claimed by Alexander and Gulyaev (between 
the vertical bars in their Fig 6) it would appear 
that the goodness of the fit exists only because 
the widths, or shapes, of the profiles generated by 
these authors are incorrect. 

3. Discussion 

The line temperature s were also examined by 
Alexander and GulvaevI (|2012l see their Fig 5). 
For RMS noise levels between 0.5 and 1 mK, which 
was the range we obtained for our telescope data, 
their line temperature vs n curve falls significantly 
below our Orion line strengths above n = 225. We 
also found that the telescope line areas decreased 
more slowl y than the test lin e areas above (n. An) 
= (202,8) (|Bell et al. II2OOOI see their Figs 3 and 



5). It was later pointed out (jBellll201l[ ) that this 
could be explained if the widths of these Orion 
lines were intrinsically narrower than predicted. 
In fact, this was considered as further proof that 
the apparent line narrowing above n = 200 was 
intrinsic and unrelated to the processing. This is 
related to the fact that in overlapping frequency- 
switched lines the amount of power lost due to 
processing decreases as the lines get narrower. 

Although these authors have found a possible 
explanation for line narrowing at high n, our previ- 
ous conclusion is not convincingly altered by their 
work since, (a) they have no independently ob- 
served telescope lines to compare their processed 
generated lines to, (b) the processed widths of the 
generated lines they are working with do not fit 
our telescope data below n = 200, nor (c) above 
200 after using their noise/no noise reduction fac- 
tors, and, d) for 1 mK noise, the temperatures of 
their processed lines cannot explain the measured 
line strengths we obtain above n ^ 220. 

These authors argue that the line narrowing we 
see in the telescope data is the result of " forcing" 
a program to fit Gaussians to features dominated 
by noise. Although noise on a line profile can in- 
crease the uncertainty in its measured width it is 
highly unlikely that this could cause such a con- 
sistent and rapid narrowing during the Gaussian 
fitting process when there are as many samples 
across the profile as here. They could have easily 
determined whether it was the Gaussian fitting in 
the presence of a decreasing s/n ratio that caused 
the narrowing by using their processed generated 
profiles without noise. By first fitting Gaussians 
to these profiles and then repeating the fitting af- 
ter noise was added to the same processed profiles, 
they could have easily determined whether it was 
the fitting in the presence of noise that caused the 
narrowing. 

4. Future Work 



Recent analysis by iHevI (|2012l) may have some 
bearing on the discre pancies betwe e n what impact 
broadening theories ( GriemI 1967 : WatsonI 2006 ) 
predict and the narrowing that has been seen in 
high-n Rydberg-Rydberg recombination spectra 
from Galactic HII regions. His results suggest that 
these theories have not allowed for the possibly sig- 
nificant role of the plasma (ion) microfield back- 
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ground in perturbing the radiator states which en- 
ter the impact broadening calculations. Violation 
of the parity selection rule would result in some 
net diminution per perturber in the effectiveness 
of the collisions of either type (electron or ion) 
in " allowed" transitions between sub-states of the 
same principal quantum number. In spite of the 
production of lin e str e ngths for 'optically forbid- 
den' transitions, iHeyl (|2012h finds that an over- 
all loss in collision strength takes place, according 
to the line strength sum rules, because of appre- 
ciable 'self-strength' contributions which would be 
absent for atomic eigenstates unperturbed by the 
microfield background. For sufficiently large prin- 
cipal quantum number, a significant departure of 
line widths from a linear dependence on electron 
(proton) density may thereby be indicated. 

Although a theoretical explanation for the ap- 
parent line n arrowing a t high n may eventually be 
forthcoming (|Hevll2012l ). the main conclusion that 
can be drawn from the results reported here is 
that attempting to solve this line-narrowing mys- 
tery using multiple-overlap frequency switching is 
unlikely to yield conclusive results. It is therefore 
suggested that using a completely different ob- 
serving technique may be the best way to address 
the problem, at least empirically (see for example 
Roshi et al., 2012). This is especially true when 
the problem is now reduced simply to one of de- 
termining whether the linewidths in question are 
stronger and narrower than predicted. Stronger, 
narrower lines are much more easily detected in 
the presence of baseline structure than are weak, 
wide lines. 

In light of this, the following is suggested as 
an approach that should give more definitive re- 
sults for rt-values above n — 200. Because the 
predicted, raw impact-broadened linewidths near 
n = 220 are so large (FWHM > 90 km/s), and 
the widths implied if the line narrowing is real are 
so much smaller (FWHM ~ 25 km/s) it should be 
easy to determine which is correct using position- 
switching together with a sufficiently long integra- 
tion time. Selecting a line, or lines, that fall in re- 
gions of the spectrum free from other strong con- 
fusing lines would be a pre-requisite. Also, posi- 
tion switching between Orion and another strong 
continuum source that does not contain recom- 
bination lines should allow most of the baseline 
structure to be cancelled. This way questions like 



how noise on the generated spectrum might affect 
the result, or whether or not the comparison line 
profiles have been generated correctly, do not en- 
ter the picture. 

5. Conclusion 

Using a frequency switching and multiple over- 
lap observing technique we have found evidence 
that above n =200 (near 6 GHz) the widths of 
impact broadened Hydrogen recombination lines 
may get much narrower and stronger than pre- 
dicted. Although Alexander and Gulyaev have 
claimed that this narrowing occurs because noise 
was present on the telescope line profiles it has 
been demonstrated here that because of errors in 
their analysis their claim may not be valid. If so, 
and the line narrowing is real, this means that 
the lines in question should be much easier to de- 
tect using more conventional observing techniques 
than previously thought. Therefore, although a 
theoretical explanation for the line narrowing may 
yet be found, for the moment the logical conclu- 
sion is that the best way to solve this mystery 
is to observe one or more of these lines using a 
different observing technique such as single-dish 
position switching, or an interferometer. 
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